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We lay out an experiment to realize time-reversal invariant topological insulators in alkali atomic gases. We 
introduce an original method to synthesize a gauge field in the near-field of an atom-chip, which effectively 
mimics the effects of spin-orbit coupling and produces quantum spin-Hall states. We also propose a feasible 
scheme to engineer sharp boundaries where the hallmark edge states are localized. Our multi-band system 
has a large parameter space exhibiting a variety of quantum phase transitions between topological and normal 
insulating phases. Due to their remarkable versatility, cold-atom systems are ideally suited to realize topological 
states of matter and drive the development of topological quantum computing. 
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Topological insulators are a broad class of unconventional 
materials that are insulating in the interior but conduct along 
the edges. The edge transport is topologically protected and 
dissipationless. This subject has emerged as a new frontier, 
discovering novelties in the single-particle band theory and 
providing a new impetus to the many -body physics of strongly 
correlated systems. Until recently, the only known topological 
insulators - quantum Hall (QH) states - violated time-reversal 
(TR) symmetry. However, the discovery of the quantum spin 
Hall (QSH) effect demonstrated the existence of novel topo- 
logical states not rooted in time-reversal violations (TJ-O, and 
has opened the possibility to design new spintronic devices ex- 
ploiting the spin-dependent currents carried by the edge states. 

Realizing topological insulators with cold atoms is partic- 
ularly attractive, and setups exploring the TR-breaking case 
have been envisaged (U 0. In this Letter, we propose a 
concrete setup using fermionic 6 Li subjected to a synthetic 
gauge field, which provides an archetypical system for in- 
vestigating the QSH phase. There are numerous proposals 
for engineering gauge fields |6], which generally depend on 
laser-induced Raman coupling between internal atomic states. 
Such a method was recently implemented (7) for bosonic 87 Rb 
atoms, but would lead to large spontaneous emission rates for 
the alkali fermions. Here we describe a setup that combines 
state-independent optical potentials with micron- scale state- 
dependent magnetic potentials in an atom chip. By completely 
eliminating spontaneous emission, this approach makes prac- 
tical the realization of gauge fields for all alkali atoms. We 
demonstrate that these synthetic fields lead to the purest real- 
ization of the QSH effect and allow to explore striking aspects 
of this topological state of matter. In particular, the stability 
of the topological phases against interactions can be explored 
by means of Feschbach resonances. 

In condensed-matter systems, the QSH effect originates 
from a material's intrinsic spin-orbit coupling (TJ 13. Note 
that such a coupling is analogous to a non-Abelian gauge 



field A, = A<r z , where <r x ,y,z are the Pauli matrices. This 
observation emphasizes that the QSH effect consists of spin- 
1/2 fermions where the two spin components are described as 
QH states at equal but opposite "magnetic fields." We hereby 
demonstrate how to synthesize such a gauge field in an optical 
lattice and show how it leads to QSH physics. 

Our proposal for realizing a fermionic model with a SU(2) 
gauge structure requires four atomic states \gi) = \F = 
1/2, m F = 1/2), \g 2 ) = |3/2,-l/2), \e±) = |3/2, 1/2), 
and \e2) = 1 1/2, — 1/2), in a square lattice described by the 
Hamiltonian 

7~L = — t ^ c L+l,n eZ6>XC ™,™ + C ^m,n+l eldy C m,n + n -C 

m,n 

+ A st ag ^]( — I)™ c ln,n C m,n' (1) 

Cm in is a 2-component field operator defined on a lattice site 
(x = ma, y = na), a is the lattice spacing, m, n are inte- 
gers, and t is the nearest-neighbor hopping. In this model, 
the Peierls phases # x ,y result from a synthetic gauge field 
that modifies the hopping along x and y. Here, all the states 
experience a primary lattice potential V\{x) = V x sm 2 (kx) 
along x which gives rise to a hopping amplitude t « 0.4 kHz. 
A secondary much weaker lattice V2(x) = 2A stag siri 2 (kx/2) 
slightly staggers the primary lattice with A sta g ~ t. These lat- 
tices, with an approximate period of a = 2 /im, are produced 
by two pairs of A = 1064 nm lasers, slightly detuned from 
each other and incident on the atom chip's reflective surface 
[Fig. [TJa)]- Additionally, these beams create a lattice along z 
with a 0.55 /im period, confining the fermions to a 2D plane. 
We study the SU(2) hopping operators 

# x = 27T7<j x , Oy — 27rxaa z , (2) 

where we set a = H — 1. The hopping operator 6 y corre- 
sponds to opposite "magnetic fluxes" ±a for each spin com- 
ponent, whereas # x mixes the spins. Our setup thus provides 
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a SU(2) generalization of the well-known two-dimensional 
electron gas in a magnetic field |9|. In order to engineer 
these state-dependent tunnelings, the states \g) and |e) must 
experience oppositely- signed lattices along y [Fig.[ljb)]. This 
can be directly implemented with the Zeeman shift g/ns \B\ 
of atoms provided that ^1,2) and \ei^) have equal, but op- 
posite magnetic moments g. In our proposal, the magnetic 
moments are correctly signed and differ by less than 1% in 
magnitude at a bias field B = 0.25 G. With these states, a 
state-dependent lattice potential can be generated by an array 
of current-carrying wires with alternating +7 and —7 currents, 
spaced by a distance a. A modest 7 = 5 jiA current (8) in wires 
3 /im below the chip-surface produces a 6 El Zeeman lattice 
[El = h 2 /8Ma is the lattice recoil energy], with a negligi- 
ble 3 Hz hopping matrix element. The assisted-hopping along 
y, with an ^-dependent phase, can be realized with an addi- 
tional grid of wires spaced by a = 2 /im along x, with currents 
7 m . This provides moving Zeeman lattices with wave- vector 
g, where I m = Io sm(qma — ut), leading to effective "Raman 
couplings". The uj/2tt w 228 ± 0.23 MHz transitions indi- 
cated with arrows in Fig. 1 (b), are independently controllable 
in phase, amplitude, and wave-vector by commanding con- 
current running waves at the indicated resonant frequencies. 
The minimum wavelength d = 2ir/q of this moving lattice is 
Nyquist limited by d > 2a. In the frame rotating at the angular 
frequency oj, and after making the rotating wave approxima- 
tion the coupling terms have the desired form t ex.p(iqma). A 
potential gradient along y detunes this Raman coupling into 
resonance and is produced by simply shifting the center of 
the harmonic potential. In the model Hamiltonian ([T])-([2]), the 
phase for hopping along y is a = qa/2n = a/ dm terms of 
physical parameters. Our scheme also requires a contribution 
to the hopping along x that mixes the |e) and \g). This can 
be realized using a Zeeman lattice moving along x, but tuned 
to drive transitions between \gi) — >> |^) and \ei) — » |e2). 
The corresponding control parameter 7, and the additional 
staggered potential which induces alternate on-site energies 
e = ±A sta g along the x-axis, are shown below to enrich the 
QSH physics in a novel manner. Note that the most sensitive 
parameter is the resonance condition required for the "Raman 
couplings". Its stability relies on the absolute control of the 
state-dependent potentials provided by the RF magnetic fields. 

The engineered Hamiltonian ([l}-([2]) satisfies TR invari- 
ance, since it commutes with the TR-operator defined as 
T = icr y K, where K is the complex-conjugate operator. This 
synthetic, yet robust, TR symmetry enables the realization of 
^-topological insulators in cold-atom laboratories. We stress 
that the two components of the field operators correspond in 
general to a pseudo-spin 1/2, but in our proposal refer specifi- 
cally to spin components of 6 Li in its electronic ground state. 

In the absence of the confining trap Konf(#> y), the system 
can be solved on an abstract cylinder. We first study the nature 
of topological insulators with this partially closed geometry 
and then show how they can be detected when the trap is ap- 
plied to the realistic open geometry. When 7 = 0, Eqs. ([l])-([2]) 
describe two uncoupled QH systems and for generic a = p/q, 




Figure 1 : (a) Experimental layout showing the origin of optical (state 
independent) and magnetic (state dependent) potentials and coupling 
fields. A state-independent, staggered, lattice along x is formed by 
the separate interference of two pairs of A = 1064 nm laser beams 
slightly detuned from each other to eliminate cross interference. The 
respective intersection angles are chosen so the lattice period differ 
by a factor of two. Both beams reflect from the chip- surface and 
form vertically aligned lattices, trapping the degenerate Fermi gas 
about 5 /im above the surface. The inset shows the chip geome- 
try, from top to bottom: a reflective chip- surface, gold wires aligned 
along y placed every 2 /im along x (producing the \g)-\e) coupling), 
and finally gold wires aligned along x with a 2 /im spacing (produc- 
ing the state dependent lattice), (b) Computed atomic potentials and 
radio-frequency (RF) driven "Raman" transitions. 

where p, q are integers, the fermion band- structure splits into 
q subbands |9|. Our setup thus provides a multi-gap system, 
where a variety of band-insulators can be reached by varying 
the atomic filling factor. As discussed below, some of these in- 
sulators are topologically non-trivial and feature gapless edge 
states. The latter are localized at the boundaries of the sample, 
and correspond to gapless excitations. When the Fermi energy 
Ef lies inside a bulk gap, the presence of these states is re- 
sponsible for the spin transport along the edges. In this Letter, 
we present results for the case a = 1/6, which exhibits ex- 
tremely rich phase diagrams with almost all possible topolog- 
ical phase transitions. The topological phases discussed be- 
low are robust against small variations 5a ~ 0.01 and rely on 
the existence of bulk gaps, which are continuously deformed 
when a is varied 0. Other configurations of the gauge field 
could be experimentally designed and would lead to similar 
effects. In Fig. [2ja), illustrating the spectrum for 7 = 0, the 
bulk bands are clearly differentiated from the gapless edge 
states within the bulk gaps. In the lowest bulk gap, the edge 
state channel A, and its TR-conjugate B, correspond to local- 
ized excitations which travel in opposite directions [Fig.[2jb)]. 
Thus the boundaries are populated by a single Kramers pair of 
counter-propagating states, each corresponding to an opposite 
spin: this lowest bulk gap describes a topological QSH phase. 
Conversely, the next gap located at E « — 1 is traversed by 
an even number of Kramers pairs: it is hence topologically 
equivalent to a normal band insulator (NI) O. 

An alternative approach to the above even-odd criteria re- 
lies on the computation of the Z2 -topological invariant v that 
characterizes the bulk gaps dEJ: KQSH) = 1 and i/(NI) =0. 
When 7 = 0, spin is conserved and the spin-conductivity is 
quantized as a s = ev/2ir (TJEl- We verified that the four 
gaps depicted in Figj2ja) are indeed associated to the sequence 
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^ = {1,0,0,1}. Surprisingly, topological insulators with dis- 
tinct orders can be realized in this multi-band scenario by sim- 
ply varying the atomic filling factor. 
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Figure 2: (a) Energy spectrum of the uncoupled system (7 = A sta g = 
and a = 1/6) computed in a cylindrical geometry: the bulk energy 
bands (thick blue bands) are traversed by edge states (thin purple 
lines), (b) Schematic representation of the two edge states, A and B, 
that lie inside the first bulk energy gap depicted in (a). The spins trav- 
eling around the edges are respectively represented by red and green 
arrows, (c) Energy bands E(k y ) for 7 = and a = 1/6, with an 
external staggered potential A sta g = 0.5£ and 1.51 The topological 
phases associated with the open bulk gaps are indicated. The pur- 
ple rectangles highlight the NI to QSH phase transition that occurs 
around Ep = —t. (d) Phase diagram in the (7, A sta g) -plane in the 
vicinity of the uncoupled case 7 = for Ep — —t. 

In our multi-band system, the lattice-potential distortions 
can drive direct transitions between normal and topological 
insulating states. Figure |5Jc) shows the bulk gaps and edge 
states for successive values of the experimentally controllable 
staggered potential A sta g. We demonstrate that A sta g induces 
a quantum phase transition (QPT) from a NI to a QSH phase 
even in the uncoupled case 7 = 0. This transition occurs within 
the bulk gaps around E = ±1 at the critical value A sta g = t. At 
half-filling, the existence of Dirac points resist the opening of 
the gap for small staggered potential and eventually lead to a 
NI phase for A sta g > 1.25t. The Z2 -index analysis provides an 
efficient tool to obtain the full phase diagram in the (7, A stag )- 
plane. The phase diagram represented in Fig. [2fd) has been 
obtained numerically by evaluating the index v inside the gap 
at E w —t, for small spin-mixing 7 < 0.1. We observe three 
distinct phases: metallic (blue), QSH (red), and NI (cyan). 
These three phases coexist at a tri-critical point situated at 7 = 
and A stag = t. The QSH phase occurs for a wide range of 7, 
indicating the robustness of this topological phase under small 
spin-mixing perturbations. 

The possibilities offered by cold-atom experiments enables 
us to consider the strong coupling regime corresponding to 



7 = 0.25. In this limit, the previously independent QH sub- 
systems (7 = 0) become maximally coupled, drastically mod- 
ifying the topological phase transitions presented above. In 
Fig. |3ja), we illustrate the bulk bands and edge states for 
Astag = 0.5£ and 1M, and the gaps are labeled according 
to the even-odd number of TR pairs. In this strong-coupled 
regime, a radically different scenario emerges: opposite phase 
transitions occur successively in the neighboring gaps. First, 
gap-closings around E w ±1 occur and trigger NI to QSH 
phase transitions at A stag = t. Then, for A stag = 1.25t, the 
opposite transition QSH-NI occurs at half-filling (E = 0). To 
fully capture the richness of this phenomenon, we numerically 
compute the index v for a wide range of the parameters around 
7 = 0.25. At half-filling, the phase diagram features tri-critical 
points, and the QSH-NI phase transitions occur along a well- 
defined curve [Fig.[3jb)]. On the other hand, in the neighbor- 
ing gaps, the NI to QSH phase transition is separated by an 
intermediate metallic region [Fig.[3jc)]. Therefore, by manip- 
ulating A stag , Bp, and the coupling 7, it is possible to explore 
different topological phase transitions within the several bulk 
gaps. These novel features endow the topological phase dia- 
gram with an intrinsic richness and complexity, not present in 
other condensed-matter realizations of the QSH effect. 
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Figure 3 : (a) Energy bands E ( k y ) for 7 = . 25 and a = 1 /6, with an 
external staggered potential A s t ag = 0.5£ and 1.5£. The topological 
phases associated to the bulk gaps are indicated, (b)-(c) Phase dia- 
grams in the (7, A sta g) -plane in the vicinity of the maximally coupled 
case 7 = 0.25 for (b) E F = and (c)E F = -t. 

We now describe a new feasible scheme to engineer a sharp 
interface where edge states can be localized. This is essential 
for detecting topological states in optical lattices, where the 
indispensable harmonic trap used to confine atoms destroys 
the edge states when V con f (edge) ~ A, A being the bulk gap's 
width [ 5 ]. The key aspect of our proposal exploits the fact that 
the hopping along the ^-direction, t y , is controlled by spatially 
periodic RF transitions and hence can be tuned. Since the har- 
monic trap has a minimal effect at the center of the trap, we 
divide the chip into three regions: the central region is char- 
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Figure 4: (a) Position dependent coupling strength. The shaded blue 
region displays the position dependent coupling strength averaged 
over one period of oscillation of the RF fields. The black trace depicts 
the Zeeman shift from the computed magnetic field along x which 
produces the RF coupling (at a representative time), (b) Discrete en- 
ergy spectrum (blue dots) and typical amplitudes y)\ 2 in the 
presence of a harmonic potential and anisotropic hopping. The open 
lattice has 42 x 42 sites, 7 = A sta g = and a = 1/6. The harmonic 
potential V CO nf(x,y) oc x 2 + y 2 is such that V CO nf(42,42) = 0.5£ 
and the hopping parameters of the inner and outer regions are respec- 
tively t' y — 2£ and = 0.1£. The amplitudes y)\ 2 mark three 
distinct regimes in the associated spectrum: the outer (light green), 
the inner (dark blue) and the edge- state (magenta) regimes. 



acterized by a large hopping t f y , while the two surrounding 
regions feature small hoppings t y <C t' y . This can be realized 
by abruptly changing the current in the "Raman" wires to a 
much smaller value, thereby reducing the coupling matrix el- 
ement on the single lattice site scale [Fig.|4fa)]. The resulting 
highly anisotropic hopping creates a sharp interface where the 
edge states of the central region - a topological insulator - 
localize. By controlling the strength of the Raman coupling, 
we squeeze the energy bands describing the outer parts so that 
they do not interfere with the bulk gaps of the central topo- 
logical phase. We show that the helical edge states are ro- 
bust and exactly localized at the designed interface. Since 
the topological phases are confined in the center of the trap, 
one verifies that the phase diagrams discussed above are valid 
for a much wider range of the harmonic potential's strength. 
Fig.Qb) illustrates the discrete energy spectrum and the typi- 
cal wave-functions marking three distinct regimes. When the 
Fermi energy lies in the edge-state regime, the edge states are 
robustly localized within the interface. These edge states are 
composed of counter-propagating spin textures, in direct ana- 
log to the spin-polarized chiral edges states of integer QH sys- 
tems. That these edge states are chiral was unambiguously 
measured by detecting a "pulse" of electrons as they propa- 
gated around - and not through - the sample ifTTTl . It is simple 



to revisit this measurement in cold atom systems by using a 
focused laser beam to launch edge excitations, followed by 
spin-sensitive imaging. To confirm the topological insulat- 
ing state, this excitation must travel along the edge with the 
two spin-components traveling in opposite directions. Further 
evidence of these edge-states could be probed through light 
Bragg scattering HUTU). When 7 = 0, the spin Hall conduc- 
tivity is given by a s = e(A^mod2)/27r, where is a Chern 
number, which can be evaluated through the Streda formula 
applied to the spin-up density lfT2l . The anisotropy t' y > £ 
leads to bulk gaps of the order A ~ £, requiring cold, though 
realistic, temperatures T~ 10 nK to detect the QSH phase. 

We have described a concrete and realistic proposal of syn- 
thetic gauge fields in optical lattices that overcomes the severe 
drawbacks affecting earlier schemes. We showed how such 
gauge fields are ideally suited to experimentally realize the 
most transparent QSH phase and hence allow to explore the 
validity of the Z 2 -classification against interaction and dis- 
order d El- In our multi-band system, a staggered poten- 
tial is shown to drive gap-dependent QPT's which constitute a 
unique and rich feature. The cold-atom realization of topolog- 
ical band-insulators and helical metals proposed in this paper 
will pave the way for engineering correlated topological su- 
perfluids and insulators. Considering atoms with more inter- 
nal states, it is possible to envisage situations where the edge 
states present a richer spin structure, and thus offer the oppor- 
tunity to explore new avenues and exotic topological phases. 
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